Bioinformatics and system biology approach to identify potential common pathogenesis for COVID-19 infection and osteoarthritis.

A growing of evidence has showed that patients with osteoarthritis (OA) had a higher coronavirus 2019 (COVID-19) infection rate and a poorer prognosis after infected it. Additionally, scientists have also discovered that COVID-19 infection might cause pathological changes in the musculoskeletal system. However, its mechanism is still not fully elucidated. This study aims to further explore the sharing pathogenesis of patients with both OA and COVID-19 infection and find candidate drugs. Gene expression profiles of OA (GSE51588) and COVID-19 (GSE147507) were obtained from the Gene Expression Omnibus (GEO) database. The common differentially expressed genes (DEGs) for both OA and COVID-19 were identified and several hub genes were extracted from them. Then gene and pathway enrichment analysis of the DEGs were performed; protein-protein interaction (PPI) network, transcription factor (TF)-gene regulatory network, TF-miRNA regulatory network and gene-disease association network were constructed based on the DEGs and hub genes. Finally, we predicted several candidate molecular drugs related to hub genes using DSigDB database. The receiver operating characteristic curve (ROC) was applied to evaluate the accuracy of hub genes in the diagnosis of both OA and COVID-19. In total, 83 overlapping DEGs were identified and selected for subsequent analyses. CXCR4, EGR2, ENO1, FASN, GATA6, HIST1H3H, HIST1H4H, HIST1H4I, HIST1H4K, MTHFD2, PDK1, TUBA4A, TUBB1 and TUBB3 were screened out as hub genes, and some showed preferable values as diagnostic markers for both OA and COVID-19. Several candidate molecular drugs, which are related to the hug genes, were identified. These sharing pathways and hub genes may provide new ideas for further mechanistic studies and guide more individual-based effective treatments for OA patients with COVID-19 infection.

Database Resources of the National Genomics Data Center, China National Center for Bioinformation in 2022.

The National Genomics Data Center (NGDC), part of the China National Center for Bioinformation (CNCB), provides a family of database resources to support global research in both academia and industry. With the explosively accumulated multi-omics data at ever-faster rates, CNCB-NGDC is constantly scaling up and updating its core database resources through big data archive, curation, integration and analysis. In the past year, efforts have been made to synthesize the growing data and knowledge, particularly in single-cell omics and precision medicine research, and a series of resources have been newly developed, updated and enhanced. Moreover, CNCB-NGDC has continued to daily update SARS-CoV-2 genome sequences, variants, haplotypes and literature. Particularly, OpenLB, an open library of bioscience, has been established by providing easy and open access to a substantial number of abstract texts from PubMed, bioRxiv and medRxiv. In addition, Database Commons is significantly updated by cataloguing a full list of global databases, and BLAST tools are newly deployed to provide online sequence search services. All these resources along with their services are publicly accessible at https://ngdc.cncb.ac.cn.

Bioinformatics Analysis of Candidate Genes and Pathways Related to Hepatocellular Carcinoma in China: A Study Based on Public Databases.

Background and Objective: Hepatocellular carcinoma (HCC) is a highly aggressive malignant tumor of the digestive system worldwide. Chronic hepatitis B virus (HBV) infection and aflatoxin exposure are predominant causes of HCC in China, whereas hepatitis C virus (HCV) infection and alcohol intake are likely the main risk factors in other countries. It is an unmet need to recognize the underlying molecular mechanisms of HCC in China. Methods: In this study, microarray datasets (GSE84005, GSE84402, GSE101685, and GSE115018) derived from Gene Expression Omnibus (GEO) database were analyzed to obtain the common differentially expressed genes (DEGs) by R software. Moreover, the gene ontology (GO) functional annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed by using Database for Annotation, Visualization and Integrated Discovery (DAVID). Furthermore, the protein-protein interaction (PPI) network was constructed, and hub genes were identified by the Search Tool for the Retrieval of Interacting Genes (STRING) and Cytoscape, respectively. The hub genes were verified using Gene Expression Profiling Interactive Analysis (GEPIA), UALCAN, and Kaplan-Meier Plotter online databases were performed on the TCGA HCC dataset. Moreover, the Human Protein Atlas (HPA) database was used to verify candidate genes' protein expression levels. Results: A total of 293 common DEGs were screened, including 103 up-regulated genes and 190 down-regulated genes. Moreover, GO analysis implied that common DEGs were mainly involved in the oxidation-reduction process, cytosol, and protein binding. KEGG pathway enrichment analysis presented that common DEGs were mainly enriched in metabolic pathways, complement and coagulation cascades, cell cycle, p53 signaling pathway, and tryptophan metabolism. In the PPI network, three subnetworks with high scores were detected using the Molecular Complex Detection (MCODE) plugin. The top 10 hub genes identified were CDK1, CCNB1, AURKA, CCNA2, KIF11, BUB1B, TOP2A, TPX2, HMMR and CDC45. The other public databases confirmed that high expression of the aforementioned genes related to poor overall survival among patients with HCC. Conclusion: This study primarily identified candidate genes and pathways involved in the underlying mechanisms of Chinese HCC, which is supposed to provide new targets for the diagnosis and treatment of HCC in China.

Classifying literature mentions of biological pathogens as experimentally studied using natural language processing.

BACKGROUND: Information pertaining to mechanisms, management and treatment of disease-causing pathogens including viruses and bacteria is readily available from research publications indexed in MEDLINE. However, identifying the literature that specifically characterises these pathogens and their properties based on experimental research, important for understanding of the molecular basis of diseases caused by these agents, requires sifting through a large number of articles to exclude incidental mentions of the pathogens, or references to pathogens in other non-experimental contexts such as public health. OBJECTIVE: In this work, we lay the foundations for the development of automatic methods for characterising mentions of pathogens in scientific literature, focusing on the task of identifying research that involves the experimental study of a pathogen in an experimental context. There are no manually annotated pathogen corpora available for this purpose, while such resources are necessary to support the development of machine learning-based models. We therefore aim to fill this gap, producing a large data set automatically from MEDLINE under some simplifying assumptions for the task definition, and using it to explore automatic methods that specifically support the detection of experimentally studied pathogen mentions in research publications. METHODS: We developed a pathogen mention characterisation literature data set -READBiomed-Pathogens- automatically using NCBI resources, which we make available. Resources such as the NCBI Taxonomy, MeSH and GenBank can be used effectively to identify relevant literature about experimentally researched pathogens, more specifically using MeSH to link to MEDLINE citations including titles and abstracts with experimentally researched pathogens. We experiment with several machine learning-based natural language processing (NLP) algorithms leveraging this data set as training data, to model the task of detecting papers that specifically describe experimental study of a pathogen. RESULTS: We show that our data set READBiomed-Pathogens can be used to explore natural language processing configurations for experimental pathogen mention characterisation. READBiomed-Pathogens includes citations related to organisms including bacteria, viruses, and a small number of toxins and other disease-causing agents. CONCLUSIONS: We studied the characterisation of experimentally studied pathogens in scientific literature, developing several natural language processing methods supported by an automatically developed data set. As a core contribution of the work, we presented a methodology to automatically construct a data set for pathogen identification using existing biomedical resources. The data set and the annotation code are made publicly available. Performance of the pathogen mention identification and characterisation algorithms were additionally evaluated on a small manually annotated data set shows that the data set that we have generated allows characterising pathogens of interest. TRIAL REGISTRATION: N/A.

Presepsin Predicts Severity and Secondary Bacterial Infection in COVID-19 by Bioinformatics Analysis.

Introduction: Novel coronavirus pneumonia (COVID-19) is an acute respiratory disease caused by the novel coronavirus SARS-CoV-2. Severe and critical illness, especially secondary bacterial infection (SBI) cases, accounts for the vast majority of COVID-19-related deaths. However, the relevant biological indicators of COVID-19 and SBI are still unclear, which significantly limits the timely diagnosis and treatment. Methods: The differentially expressed genes (DEGs) between severe COVID-19 patients with SBI and without SBI were screened through the analysis of GSE168017 and GSE168018 datasets. By performing Gene Ontology (GO) enrichment analysis for significant DEGs, significant biological processes, cellular components, and molecular functions were selected. To understand the high-level functions and utilities of the biological system, the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis was performed. By analyzing protein-protein interaction (PPI) and key subnetworks, the core DEGs were found. Results: 85 DEGs were upregulated, and 436 DEGs were downregulated. The CD14 expression was significantly increased in the SBI group of severe COVID-19 patients (P < 0.01). The area under the curve (AUC) of CD14 in the SBI group in severe COVID-19 patients was 0.9429. The presepsin expression was significantly higher in moderate to severe COVID-19 patients (P < 0.05). Presepsin has a diagnostic value for moderate to severe COVID-19 with the AUC of 0.9732. The presepsin expression of COVID-19 patients in the nonsurvivors was significantly higher than that in the survivors (P < 0.05). Conclusion: Presepsin predicts severity and SBI in COVID-19 and may be associated with prognosis in COVID-19.

Girl-Boy Twins with Developmental Delay from 16p11.2 Triplication due to Biparental Inheritance from Two Parents with 16p11.2 Duplication.

The 16p11.2 duplication is a well-known cause of developmental delay and autism, but there are only 2 previously reported cases of 16p11.2 triplication. Both of the previously reported cases exhibited tandem triplication on a 16p11.2 duplication inherited from 1 parent. We report fraternal twins presenting with developmental delay and 16p11.2 triplication resulting from inheritance of a 16p11.2 duplicated homolog from each parent. This report also reviews the overlapping features in previously published cases of 16p11.2 triplication, and possible implications are discussed.

Effective deep learning approaches for predicting COVID-19 outcomes from chest computed tomography volumes.

The rapid evolution of the novel coronavirus disease (COVID-19) pandemic has resulted in an urgent need for effective clinical tools to reduce transmission and manage severe illness. Numerous teams are quickly developing artificial intelligence approaches to these problems, including using deep learning to predict COVID-19 diagnosis and prognosis from chest computed tomography (CT) imaging data. In this work, we assess the value of aggregated chest CT data for COVID-19 prognosis compared to clinical metadata alone. We develop a novel patient-level algorithm to aggregate the chest CT volume into a 2D representation that can be easily integrated with clinical metadata to distinguish COVID-19 pneumonia from chest CT volumes from healthy participants and participants with other viral pneumonia. Furthermore, we present a multitask model for joint segmentation of different classes of pulmonary lesions present in COVID-19 infected lungs that can outperform individual segmentation models for each task. We directly compare this multitask segmentation approach to combining feature-agnostic volumetric CT classification feature maps with clinical metadata for predicting mortality. We show that the combination of features derived from the chest CT volumes improve the AUC performance to 0.80 from the 0.52 obtained by using patients' clinical data alone. These approaches enable the automated extraction of clinically relevant features from chest CT volumes for risk stratification of COVID-19 patients.

Genetics Matters: Voyaging from the Past into the Future of Humanity and Sustainability.

The understanding of how genetic information may be inherited through generations was established by Gregor Mendel in the 1860s when he developed the fundamental principles of inheritance. The science of genetics, however, began to flourish only during the mid-1940s when DNA was identified as the carrier of genetic information. The world has since then witnessed rapid development of genetic technologies, with the latest being genome-editing tools, which have revolutionized fields from medicine to agriculture. This review walks through the historical timeline of genetics research and deliberates how this discipline might furnish a sustainable future for humanity.

dbGSRV: A manually curated database of genetic susceptibility to respiratory virus.

Human genetics has been proposed to play an essential role in inter-individual differences in respiratory virus infection occurrence and outcomes. To systematically understand human genetic contributions to respiratory virus infection, we developed the database dbGSRV, a manually curated database that integrated the host genetic susceptibility and severity studies of respiratory viruses scattered over literatures in PubMed. At present, dbGSRV contains 1932 records of genetic association studies relating 1010 unique variants and seven respiratory viruses, manually curated from 168 published articles. Users can access the records by quick searching, batch searching, advanced searching and browsing. Reference information, infection status, population information, mutation information and disease relationship are provided for each record, as well as hyperlinks to public databases in convenient of users accessing more information. In addition, a visual overview of the topological network relationship between respiratory viruses and associated genes is provided. Therefore, dbGSRV offers a convenient resource for researchers to browse and retrieve genetic associations with respiratory viruses, which may inspire future studies and provide new insights in our understanding and treatment of respiratory virus infection. Database URL: http://www.ehbio.com/dbGSRV/front/.

The Sequence Read Archive: a decade more of explosive growth.

The Sequence Read Archive (SRA, https://www.ncbi.nlm.nih.gov/sra/) stores raw sequencing data and alignment information to enhance reproducibility and facilitate new discoveries through data analysis. Here we note changes in storage designed to increase access and highlight analyses that augment metadata with taxonomic insight to help users select data. In addition, we present three unanticipated applications of taxonomic analysis.

Getting Started with the IDG KMC Datasets and Tools.

The Illuminating the Druggable Genome (IDG) consortium is a National Institutes of Health (NIH) Common Fund program designed to enhance our knowledge of under-studied proteins, more specifically, proteins unannotated within the three most commonly drug-targeted protein families: G-protein coupled receptors, ion channels, and protein kinases. Since 2014, the IDG Knowledge Management Center (IDG-KMC) has generated several open-access datasets and resources that jointly serve as a highly translational machine-learning-ready knowledgebase focused on human protein-coding genes and their products. The goal of the IDG-KMC is to develop comprehensive integrated knowledge for the druggable genome to illuminate the uncharacterized or poorly annotated portion of the druggable genome. The tools derived from the IDG-KMC provide either user-friendly visualizations or ways to impute the knowledge about potential targets using machine learning strategies. In the following protocols, we describe how to use each web-based tool to accelerate illumination in under-studied proteins. © 2022 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Interacting with the Pharos user interface Basic Protocol 2: Accessing the data in Harmonizome Basic Protocol 3: The ARCHS4 resource Basic Protocol 4: Making predictions about gene function with PrismExp Basic Protocol 5: Using Geneshot to illuminate knowledge about under-studied targets Basic Protocol 6: Exploring under-studied targets with TIN-X Basic Protocol 7: Interacting with the DrugCentral user interface Basic Protocol 8: Estimating Anti-SARS-CoV-2 activities with DrugCentral REDIAL-2020 Basic Protocol 9: Drug Set Enrichment Analysis using Drugmonizome Basic Protocol 10: The Drugmonizome-ML Appyter Basic Protocol 11: The Harmonizome-ML Appyter Basic Protocol 12: GWAS target illumination with TIGA Basic Protocol 13: Prioritizing kinases for lists of proteins and phosphoproteins with KEA3 Basic Protocol 14: Converting PubMed searches to drug sets with the DrugShot Appyter.

Spatio-temporal dynamics of intra-host variability in SARS-CoV-2 genomes.

During the course of the COVID-19 pandemic, large-scale genome sequencing of SARS-CoV-2 has been useful in tracking its spread and in identifying variants of concern (VOC). Viral and host factors could contribute to variability within a host that can be captured in next-generation sequencing reads as intra-host single nucleotide variations (iSNVs). Analysing 1347 samples collected till June 2020, we recorded 16 410 iSNV sites throughout the SARS-CoV-2 genome. We found ∼42% of the iSNV sites to be reported as SNVs by 30 September 2020 in consensus sequences submitted to GISAID, which increased to ∼80% by 30th June 2021. Following this, analysis of another set of 1774 samples sequenced in India between November 2020 and May 2021 revealed that majority of the Delta (B.1.617.2) and Kappa (B.1.617.1) lineage-defining variations appeared as iSNVs before getting fixed in the population. Besides, mutations in RdRp as well as RNA-editing by APOBEC and ADAR deaminases seem to contribute to the differential prevalence of iSNVs in hosts. We also observe hyper-variability at functionally critical residues in Spike protein that could alter the antigenicity and may contribute to immune escape. Thus, tracking and functional annotation of iSNVs in ongoing genome surveillance programs could be important for early identification of potential variants of concern and actionable interventions.

Identification of hub genes associated with COVID-19 and idiopathic pulmonary fibrosis by integrated bioinformatics analysis.

INTRODUCTION: The coronavirus disease 2019 (COVID-19), emerged in late 2019, was caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The risk factors for idiopathic pulmonary fibrosis (IPF) and COVID-19 are reported to be common. This study aimed to determine the potential role of differentially expressed genes (DEGs) common in IPF and COVID-19. MATERIALS AND METHODS: Based on GEO database, we obtained DEGs from one SARS-CoV-2 dataset and five IPF datasets. A series of enrichment analysis were performed to identify the function of upregulated and downregulated DEGs, respectively. Two plugins in Cytoscape, Cytohubba and MCODE, were utilized to identify hub genes after a protein-protein interaction (PPI) network. Finally, candidate drugs were predicted to target the upregulated DEGs. RESULTS: A total of 188 DEGs were found between COVID-19 and IPF, out of which 117 were upregulated and 71 were downregulated. The upregulated DEGs were involved in cytokine function, while downregulated DEGs were associated with extracellular matrix disassembly. Twenty-two hub genes were upregulated in COVID-19 and IPF, for which 155 candidate drugs were predicted (adj.P.value < 0.01). CONCLUSION: Identifying the hub genes aberrantly regulated in both COVID-19 and IPF may enable development of molecules, encoded by those genes, as therapeutic targets for preventing IPF progression and SARS-CoV-2 infections.

Coronavirus GenBrowser for monitoring the transmission and evolution of SARS-CoV-2.

Genomic epidemiology is important to study the COVID-19 pandemic, and more than two million severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genomic sequences were deposited into public databases. However, the exponential increase of sequences invokes unprecedented bioinformatic challenges. Here, we present the Coronavirus GenBrowser (CGB) based on a highly efficient analysis framework and a node-picking rendering strategy. In total, 1,002,739 high-quality genomic sequences with the transmission-related metadata were analyzed and visualized. The size of the core data file is only 12.20 MB, highly efficient for clean data sharing. Quick visualization modules and rich interactive operations are provided to explore the annotated SARS-CoV-2 evolutionary tree. CGB binary nomenclature is proposed to name each internal lineage. The pre-analyzed data can be filtered out according to the user-defined criteria to explore the transmission of SARS-CoV-2. Different evolutionary analyses can also be easily performed, such as the detection of accelerated evolution and ongoing positive selection. Moreover, the 75 genomic spots conserved in SARS-CoV-2 but non-conserved in other coronaviruses were identified, which may indicate the functional elements specifically important for SARS-CoV-2. The CGB was written in Java and JavaScript. It not only enables users who have no programming skills to analyze millions of genomic sequences, but also offers a panoramic vision of the transmission and evolution of SARS-CoV-2.

A comprehensive overview of identified mutations in SARS CoV-2 spike glycoprotein among Iranian patients.

Since late 2019, when SARS-CoV-2 was reported at Wuhan, several sequence analyses have been performed and SARS-CoV-2 genome sequences have been submitted in various databases. Moreover, the impact of these variants on infectivity and response to neutralizing antibodies has been assessed. In the present study, we retrieved a total number of 176 complete and high-quality S glycoprotein sequences of Iranian SARS-COV-2 in public database of the GISAID and GenBank from April 2020 up to May 2021. Then, we identified the number of variables, singleton and parsimony informative sites at both gene and protein levels and discussed the possible functional consequences of important mutations on the infectivity and response to neutralizing antibodies. Phylogenetic tree was constructed to represent the relationship between Iranian SARS-COV2 and variants of concern (VOC), variants of interest (VOI) and reference sequence. We found that the four current VOCs - Alpha, Beta, Gamma and Delta - are circulated in different regions in Iran. The Delta variant is notably more transmissible than other variants, and is expected to become a dominant variant. However, some of the Delta variants in Iran carry an additional mutation, namely E1202Q in the HR2 subdomain that might confer an advantage to viral/cell membrane fusion process. We also observed some more common mutations such as an N-terminal domain (NTD) deletion at position I210 and P863H in fusion peptide-heptad repeat 1 span region in Iranian SARS-COV-2. The reported mutations in the current project have practical significance in prediction of disease spread as well as design of vaccines and drugs.

Immunoinformatic Analysis Reveals Antigenic Heterogeneity of Epstein-Barr Virus Is Immune-Driven.

Whole genome sequencing of Epstein-Barr virus (EBV) isolates from around the world has uncovered pervasive strain heterogeneity, but the forces driving strain diversification and the impact on immune recognition remained largely unknown. Using a data mining approach, we analyzed more than 300 T-cell epitopes in 168 published EBV strains. Polymorphisms were detected in approximately 65% of all CD8+ and 80% of all CD4+ T-cell epitopes and these numbers further increased when epitope flanking regions were included. Polymorphisms in CD8+ T-cell epitopes often involved MHC anchor residues and resulted in changes of the amino acid subgroup, suggesting that only a limited number of conserved T-cell epitopes may represent generic target antigens against different viral strains. Although considered the prototypic EBV strain, the rather low degree of overlap with most other viral strains implied that B95.8 may not represent the ideal reference strain for T-cell epitopes. Instead, a combinatorial library of consensus epitopes may provide better targets for diagnostic and therapeutic purposes when the infecting strain is unknown. Polymorphisms were significantly enriched in epitope versus non-epitope protein sequences, implicating immune selection in driving strain diversification. Remarkably, CD4+ T-cell epitopes in EBNA2, EBNA-LP, and the EBNA3 family appeared to be under negative selection pressure, hinting towards a beneficial role of immune responses against these latency type III antigens in virus biology. These findings validate this immunoinformatics approach for providing novel insight into immune targets and the intricate relationship of host defense and virus evolution that may also pertain to other pathogens.

Targeting genomic SARS-CoV-2 RNA with siRNAs allows efficient inhibition of viral replication and spread.

A promising approach to tackle the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) could be small interfering (si)RNAs. So far it is unclear, which viral replication steps can be efficiently inhibited with siRNAs. Here, we report that siRNAs can target genomic RNA (gRNA) of SARS-CoV-2 after cell entry, and thereby terminate replication before start of transcription and prevent virus-induced cell death. Coronaviruses replicate via negative sense RNA intermediates using a unique discontinuous transcription process. As a result, each viral RNA contains identical sequences at the 5' and 3' end. Surprisingly, siRNAs were not active against intermediate negative sense transcripts. Targeting common sequences shared by all viral transcripts allowed simultaneous suppression of gRNA and subgenomic (sg)RNAs by a single siRNA. The most effective suppression of viral replication and spread, however, was achieved by siRNAs that targeted open reading frame 1 (ORF1) which only exists in gRNA. In contrast, siRNAs that targeted the common regions of transcripts were outcompeted by the highly abundant sgRNAs leading to an impaired antiviral efficacy. Verifying the translational relevance of these findings, we show that a chemically modified siRNA that targets a highly conserved region of ORF1, inhibited SARS-CoV-2 replication ex vivo in explants of the human lung. Our work encourages the development of siRNA-based therapies for COVID-19 and suggests that early therapy start, or prophylactic application, together with specifically targeting gRNA, might be key for high antiviral efficacy.